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Semiclassical transformation theory implies an integral representation for stationary-state wave 
functions ipmiq) in terms of angle-action variables {9, J). It is a particular solution of Schrodinger's 
time- independent equation when terms of order and higher are omitted, but the pre-exponential 
factor A{q, 6) in the integrand of this integral representation does not possess the correct dependence 
on q. The origin of the problem is identified: the standard unitarity condition invoked in semiclassical 
transformation theory does not fix adequately in A[q, 9) a factor which is a function of the action J 
written in terms of q and 9. A prescription for an improved choice of this factor, based on succesfuUy 
. . . reproducing the leading behaviour of wave functions in the vicinity of potential minima, is outlined, 

xj , Exact evaluation of the modified integral representation via the Residue Theorem is possible. It 

■ yields wave functions which are not, in general, orthogonal. However, closed-form results obtained 
' after Gram-Schmidt orthogonalization bear a striking resemblance to the exact analytical expressions 

(N ■ 

for the stationary-state wave functions of the various potential models considered (namely, a Poschl- 
5— ( . Teller oscillator and the Morse oscillator). 

<\ . . 

Semiclassical transformation theory approximates various quantum probability amplitudes in terms of the gener- 

■ ating functions of classical canonical transformations 0] . The theory is aesthetically very satisfying in the manner in 
' ' which it exploits parallels between quantum mechanics and classical mechanics. It has also proved eminently useful. 
I , The basic relations for amplitudes in terms of generating functions, although more than 30 years old, are still central 

^ ' to current research 

p\| , A formal application of semiclassical transformation theory implies that the bound states of a conservative system 

• of one degree of freedom have actions Jm {"m — 0,1,2,.. .) given by the appropriate Bohr-Sommcrfeld quantization 

, condition and approximate configuration space wave functions 

O . f 

^ , i^miq) = Nm / A{q,e)exp[{i/h)F,i{q,e)](j,,n{0)de, (1) 

o : ^ 

where Fci{q,9) is the generating function of the first kind (see chapter 9 in 0|) for the canonical transformation from 
Oh, conjugate cartesian variables {q,p) to the angle-action variables {9, J) for the system, the pre-exponential factor 



C^; A{q,9)=Au{q,0)^ 



>< 



-1 d^Fr, 



2T:ih dqd9 



-,1/2 

(2) 



(j)m{9) = {2T:ih)^^/^ exp{iJm9/ti) and Nm is a normalization constant j^. Equations |^ and ((SJ reduce to an 
acceptable result (namely, the WKB approximation) if the integral is evaluated in the h ^ stationary phase 
approximation. However, if one attempts to go beyond the stationary phase approximation as contemporary studies 
' do, then there are problems. As specialization of Eqs. Q and Q to various examples reveals (see below and Appendix 
C.2 in J(|), the q dependence of Au is typically spurious. In particular, wave functions do not have well-defined parity 
under circumstances when this symmetry is expected. 

The integral representation of Eq. ^ with the pre-exponential factor in Eq. |(2Jl is a particular solution of 
Schrodinger's time-independent equation when terms of order equal to or greater than ti^ are discarded, but the 
general solution (to this order) only constrains the pre-exponential factor A{q, 9) to be of the form (3| 

AG{q,9) = f{J{q,9))Au{q,9), (3) 

where / is an arbitary differcntiable function of J{q,9) — —dFci/d9 Q. I claim that, in the case of potentials 
displaying a single minimum, the factor / {J{q, 9)) can be chosen so that the corresponding approximate wave functions 
have the desired q behaviour in the vicinity of the potential's minimum. In the case of even potentials V{q) (with 
a single minimum), this is enough to guarantee that wave functions have the right parity. No matter what choice 
of / is made, the pre-exponential factors Au{q,9) and ^^(g,^) are indistinguishable if the ?i — > stationary phase 
approximation is invoked [because J{q,9) is then replaced by the constant Jm]- 
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The choice of pre-exponential factor Au{q,d) guarantees that, in the hmit ?i — > 0, the corresponding kernel 
f^uilj(^) = -^[/('Z: exp[(«/?i)Fc;(g, 6*)] in Eq. is an element of a unitary transformation matrix (or a quantum 
probability amplitude). Although this is an appealing property, it is by no means obvious that it is appropriate as 
there is no rigorous unitary quantum mechanical counterpart to transformations from cartesian phase space variables 
to (conventional) angle- action variables. In fact, one implication of the previous paragraph is that it is permissable 
to insist that )C{q, 0) = A{q, 9) exp[(i/?i)i^c;(<Z, 9)] is an element of a unitary transformation matrix in the limit ?i — > 0. 
The upshot of the considerations below is that this unitarity condition must be supplemented by more information 
to fix the pre-exponential factor completely. 

Support for the assertion above about the role of the factor f{J{q, 0)) comes from consideration of some analytically 
soluble models: the simple harmonic oscillator (discussed in Q), a Poschl- Teller oscillator described by the potential 
V{q) = Vbtan2(7rg/a) {Vq > 0) and the Morse oscillator with potential V{q) = D[l - exp(-g/d)]2 {D,d > 0). The 
Poschl- Teller and Morse models, which are taken up in the present work, also force one to consider issues related to 
the contour of integration in Eq. 

For the example of the harmonic oscillator, eigenfunctions can be obtained by integrating over the full range of 
(modulo 27r) for which Fci{q, 9) exists for the value of q under consideration Despite the apparent reasonableness of 
this prescription, it is inappropriate in the case of the Poschl- Teller potential. The corresponding limits of integration 
would contain a q-dependence |Toj incompatible with the requirement that the resulting integral representation for 
eigenfunctions satisfies Schrodinger's time-independent equation when terms of order and higher are dropped. (The 
q-dependent limits would give rise to error terms of order Ti.) 

What alternatives are there? The general formal considerations of Q suggest that the kernel K,{q,6) — 
A{q,0) exp[{i/h)Fci{q.,0)] in Eq. can be thought as a generating function for (approximate) wave functions in 
the sense that they are determined (up to multiplicative constants) by the coefficients of an expansion of )C{q,9) in 
powers of w = e~*^. More precisely, if the Bohr-Sommerfeld quantization condition Jm = (m -|- iJ,/4:)h applies {m is a 
non-negative integer and fi is the Maslov index), then, according to Eq. (10) in 0, the expansion reads 



IC{q,9{w)) ^w^^l^ 



m— 



-W 



(4) 



where 0[yj) = ilnw [with \'w\ chosen so that the series in Eq. I^J converges] and the n^'s are constants [introduced 
in the denominator to simplify Eqs. (O and ©]. The expansion in Eq. Q implies that approximate wave functions 
can be written as the Schlafli-like contour integral 



m+l+^,./4 



A {q, e{w)) exp [{i/h)F,i {q, e{w))] dw, 



(5) 



where the (counter-clockwise) contour of integration is a simple closed curve around a pole of order to -I- 1 at w = 
which does not enclose any other singularities of the integrand. 

Despite appearances, the integral representations in Eqs. 1^ and lO are closely related: since Jm = (to -f /i/4)?i in 
4>m(0), the integrand of Eq. ^ transforms into the integrand of Eq. (|3J) under the replacement of 6* by w = e^'^ as the 
variable of integration. Thus, implicit in Eq. Q is a choice of the contour of integration for the integral representation 
in Eq. If possible, it should pass through the point(s) of stationary phase which yield the WKB approximation. 
I shall not attempt to be more explicit here, because, for the purposes of this paper, the Rodrigues' formula implied 
by Eq. Q [or the Residue Theorem applied to Eq. ©J, namely 



^m(g) 



d 



7n 



n. 

to! dw 



w^^'/^K.{q,e{w)) 



(6) 



w=0 



is more useful than either of the above integral representations. 

When one can distinguish between states of different parity (as in the Poschl- Teller model), it is necessary to 
introduce separate kernels Kg for the positive {g = Q) and negative parity [g — 1) states [obtained by making different 
choices of the factor / (J(g,0))]. An expansion in powers of z = e~*^* is now appropriate, the expansion being 



JCg {q, 0{z)) = V t^^ilM, 

m— ^ 



(7) 



The associated Rodrigues' formula expresses tp2m+giq) in terms of the mth [not the (2m + g)t]i] partial derivative of 
2;-(A'/4+e)/2/(;^ (q, 0{z)) with respect to z (evaluated a.t z — 0). A division into two classes of states can also be helpful 
for models which do not respect parity (see the discussion of the Morse oscillator below). 
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The harmonic osciUator problem furnishes evidence that expansions hke that in Eq. Q or Eq. Q do, in fact, exist. 
Consistent with Eq. (20) in the kernels of definite parity are 

JC,{q,9)=yB{cos9r^^/^+e^exp(~'-yHan9] , (8) 



where the dimensionless variable y — ^ muj/fi q (m is the mass of the oscillator and lu its angular frequency) . The 
expression which results on the substitution of e~*^^ by z can be recast into the form 



IC, = 2V2+e^i/4+,/2g-,V2 ( ) y^^e^y" (9) 



suitable for expansion in powers of z Using the Maclaurin series in z for y^e^^^ and the fact that the Hermite 

1 aL 

4 dy^ 



polynomial Hn{y) = 2" exp(— -H-a)?/" (see Eq. (1) in 12]), Eq. © implies that 



where the iy9„(g)'s are the exact normalized energy eigenfunctions of the simple harmonic oscillator (iy9„ — 
[ma;/(7r?i)]^/'*(2"n!)^/^iJ„(j/)e~^ 1"^). Since the Maslov index ^ = 2 for the harmonic oscillator, Eq. ifTIHl constitutes 
an explicit realization of Eq. (O . 

The Poschl- Teller oscillator is instructive, because unlike the example of the harmonic oscillator, the results obtained 
are not exact and so it is possible to assess the quality of the approximation based on Eqs. 1^ and (jj)). Although 
the relevant transformation from cartesian phase variables to angle-action variables appears in textbooks (see, for 
example, chapter 7 in [T^l. the corresponding generating function of the first kind is not so readily available. The 
calculation of this generating function as the Legendre transform of Hamilton's characteristic function T4^(g, J) for 
the Poschl- Teller problem is straightforward but lengthy, it being necessary to treat separately each branch of the 
multiple-valued Wijq^ J) for both signs of q. Fortunately, the result can be compactly expressed [in a form suitable 
for use in the Rodrigues' formula implied by Eq. ((JJ)] as follows: 



Fciiq,0)^iJs\rL- 



-^29y + 4cos2(7rq/a)e-'2e + (i - e"^*^) 



[2 cos(7rg/a)] ^ (11) 



where the parameter Jg = av^2mVo^/7r sets the scale for classical actions (m is again the mass of the oscillator). 
In Eq. Hll|l . which holds for both signs of q and the pertinent ranges of 9 (see [l0j)i it is understood that the 
positive branch of the square root and the principal value of the logarithm function are to be adopted [l4j . Partial 



prii ] 

differentiation of Fci{q,0) with respect to 9 [13 yields 

J{q,9) ^ = J, |(l - [sin(7rq/a)/cos0]')"'^' - l| , (12) 

implying that, for the present problem, the factor f{J{q,9)) in Ag is, in effect, a function of the combination 
sin(7rq/a)/ cos9. Likewise, the second partial derivative of Fd required in Aq (and Ajj) 

= J.cos Wa — \ 1 - sin Wa /cos6if \ 13 

a cos^ 9 1 J 

amounts to the product of cot(7rg/a) with a function of sin(7rq/a)/ cos0. Pre-exponential factors must, therefore, 
be of the form [cot(7rg/a)]^/2G(sin(7r(7/a)/ cos0), where G is some differentiable function. This functional form is 
compatible with choices which possess the leading small q behaviour expected of energy eigenfunctions of definite 
parity, the simplest being 

A,iq, 9) = [cot(.,/a)] V2 ^ [eos(.,/«)] V2 i^^^!!^, (14) 

\ cos 6^ / (cos6')^'^+£' 

where, as above, g = {1) for positive (negative) parity states. Equation 1)13(1 serves to illustrate that the pre- 
exponential factor Ajj would have the wrong leading behaviour for small q. 
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With the pre-exponential factors of Eq. H14|l . the Rodrigues' formula for approximate energy eigenfunctions implies 
that 



■02m+e(?) oc l^cos (sin(7rq/a)) 



where K = Jg/Ti and 



dz'' 



(1 + _ ^ + ^(1 + z)2 - 4x2z) " 



(15) 



(16) 



z=0 



is a polynomial of order 2k + g and well-defined parity [Pj.^\—x) — {—1)^pI.^\x)]. Calculation shows that these 
approximate wave functions are not orthogonal although the degree of non-orthogonality is slight if A ^ 1 . In 
terms of x = sm^irq/a), orthogonalization amounts to the construction of a set of orthogonal polynomials on the 

interval — 1 < a: < -1-1 with weight function (1 — x^)^. This is the set of ultraspherical polynomials C^^^^^\x) 
Accordingly, after orthonormalization, the approximate wave functions are 



iV„(A) [cos{nq/a)f+'^^ Ci^+^'^^ (sin(^q/a)) 



(17) 



where the iV„(A)'s are normalization constants (n = 0, 1, 2, . . .). Remarkably, the exact energy eigenfunctions are also 
given by Eq. (^Tj) with A replaced by ^A2 + 1/4. 

The analytic expressions for the non-orthogonal wave functions in Eq. 11511 and the orthogonalized wave functions 
in Eq. (|17|l do not appear to be significantly different: barring the overall normalization factors, the coefficients of 
each power of a; = svii{'Kq/a) agree to leading order in A (see [Tg for some examples). Nevertheless, orthogonalization 
is crucial. Numerical studies reveal that it improves substantially the agreement between approximate and exact 
eigenfunctions (so that it is good provided A^ 3> 1). Perhaps more telling is the fact that, even when A ^ 1, the 
non-orthogonal wave functions can be inferior to (approximately) normalized WKB wavefunctions 19] inside the 
classically- allowed region (well away from the classical turning points). 

Analysis of the Morse oscillator serves to illustrate that it is possible to treat models which do not possess wave 
functions of good parity. The guiding principle I have found useful is to insist that, in the immediate vicinity 
of the potential minimum (at q = 0), the approximate Morse oscillator wave functions i^nio) resemble those of a 
harmonic oscillator of the same mass m and the appropriate angular frequency, namely lo = Jsl(md^) where the 
action Jg = \/2mDd [23|. Not only is the pre-exponential factor A{q, 9) at issue but also the analytic continuation to 
complex values of of the generating function Fci{q, 9) to be used in Eqs. (O or © and even the choice of Fdiq, 6) 
for real values of 9 |21jj. Adoption of the q > result for Fci{q, 9) is indicated, which, for < < tt, reads 



F,i{q,e) = jAJy{q,e)-{l-e-i/dY+kvccos (l-e-«/'*)/A(9,0) - 



(18) 

where A(g, 9) is the positive root of the relation q/d — ln[(l + A cos 6')/(l — A^)] and Arccos denotes the principal value 
of the inverse cosine function. It is understood that the principal values of all square roots present either explicitly 
or implicitly [in A(0, g)] are to be used. The analytic continuation (of interest) of Fci{q,9) to complex 9 is such that 
the alternative branch of the square root y/X^iq, 9) — {1 — e-il'^Y in Eq. ifTI^ is adopted (22] . 

In view of the fact that I compare Morse oscillator wave functions with those of the harmonic oscillator, it is still 
useful to divide them into a g = class {'02m('?)} and a p = 1 class {'02m-i-i('?)}, which correspond (near g = 0) to the 
positive and negative parity harmonic oscillator wave functions, respectively (m = 0, 1, . . .). Suitable pre-exponential 
factors for the two classes of wave functions are [e^^^ cos^ ^?+4(l — e^'^)]^^/"'Gg(A(g, 9)), where the simplest appropriate 



choice of Gq(x) is Gq(x) 



For this choice, the approximate energy eigenfunctions are found to be of the form 



V'„(g) cx (e-«/'*)^-"-i/2gxp(_Ae-«/'')P„(e-?/'') 



(19) 



where Pn{x) is a polynomial of order n and A= Jg/h. As in the case of the Poschl- Teller oscillator, the wave functions 
ipnil) in Eq- (|19|l are weakly non-orthogonal. After orthonormalization, the exact bound state wave functions are 
obtained. 

In summary, the examples in this paper illustrate just how much can be learnt about wave functions of simple systems 
using a modification of the approximate integral representation in Eq. Q ■ One change seeks to correct for obvious 
deficiencies of the pre-exponential factor implied by semiclassical transformation theory [namely, Au in Eq. Q]. 
Other changes may appear less welcome but are equally essential: the use of a complex-valued contour of integration 
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and, sometimes, an appropriate analytic continuation of the generating function. Despite the lack of intuitive appeal 
inherent in the use of a complex-valued contour, it does facilitate analytic evaluation of the integral representation 
(since the Residue Theorem may be invoked) . Typically, semiclassical integral representations improving on Eq. 
entail integration over classical trajectories 123. |24L l25l| . This is one obvious difference between these representations 
and the one explored here. There is another. Corrections to the pre-exponential factor have been guided in this 
work by the structure of the general solution in Eq. Q . This approach is inappropriate for integral representations 
involving integration over classical trajectories: the factor /( J(q, 0)) in Eq. (|2Jl becomes an invisible constant [since 
J{q, 9) is equal to the fixed action on a trajectory]. 



I would like to thank Ghassan I. Ghandour for bringing 12] to my attention and for suggesting investigation of the 
Poschl- Teller model. 
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